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ABSTRACT 

A three-dimensional finite element model of the chemical vapor deposition 
(CVD) of boron on a carbon substrate (B/C) is developed. The model includes 
an expansion of the boron after deposition due to atomic rearrangement and in- 
cludes creep of the boron and carbon. Curves are presented to show how the 
principal residual stresses and the filament elongation vary as the parameters 
defining deposition strain and creep are varied. The calculated results are 
compared with experimental axial residual stress and elongation measurements 
made on B/C filaments. This comparison requires that for good agreement be- 
tween calculated and experimental results, the deposited boron must continue 
to expand after deposition, and that the build-up of residual stresses is 
limited by significant boron and carbon creep rates. 

INTRODUCTION 

Boron filaments made by the chemical vapor deposition (CVD) of boron onto 
a tungsten or carbon substrate exhibit high moduli and strengths; however, a 
detailed understanding of the elongation of the filament during deposition and 
the large residual stresses in the resulting filament is yet needed. For 
boron deposited on a carbon substrate (B/C), this elongation is as much as 
5 percent (Ref. 1). Measurements of the axial residual stresses in 102 ym 
(4.0 mil) diameter B/C filament show these stresses to be large, 790 MPa 
(115 ksi) compressive stress near the surface and 290 MPa (42 ksi) tensile 
stress near the carbon core (Ref. 2). 

Witucki (Ref. 3) estimated the residual stresses by assuming that the 
boron is deposited in a stressed state and that the boron behaved elastically. 
His calculated axial residual stress distributior , although agreeing with ex- 
perimental values at the surface, increased rapidly for decreasing radial 
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location inside the filament to values much larger than are observed experi- 
mentally. Mehalso (Ref. 1) suggested that the residual stresses and the 
elongation were produced by an expansion of the boron after deposition due to 
atomic rearrangement. Eason (Ref. 4) has measured the elongation during the 
chamical vapor deposition of boron on tungsten and boron on carbon. Eason, 
using a one-dimensional model in which he assumes, as did Mehalso, that the 
boron continues to expand after deposition and that the residual stresses are 
limited to a given maximum value, calculate- 1 the filament elongation during 
boron deposition. Although he obtained ad agreement between his calculated 
and measured elongation results, his model did not give good agreement with 
the measured axial residual stresses. Because his model was one-dimensional, 
he was not able to calculate the radial and circumferential stresses. 

The objective of the work reported here was to develop a three-dimensional 
model of chemical vapor deposition of boron and to use this to calculate resid- 
ual stresses and elongation of the filament during deposition. Because the 
model developed was three-dimensional, the principal stresses in the radial, 
circumferential, and axial directions could all be calculated. In this model, 
the boron was allowed to expand after deposition as suggested by Mehalso and 
used by Eason in his model. In addition, the boron was allowed to creep ac- 
cording to a stress-dependent creep law. These two additions in this model 
(the use of three dimensions and the stress-dependent creep) produce calculated 
results for the residual stresses and filament elongation in excellent agree- 
ment with experimental measurements. 

DEPOSITION MODEL 

It is assumed in this model that the filament is always axisymmetric. 

The elements in the finite element analysis used were taken to be cylindrical 
shells as shown In Figure 1. Here, r c is the core (substrate) radius, and 
rji is the outer radius of the i^ shell or element. The nodes of the ith 
element are then located at r i-i and Tj. Because of this symmetry, it is 
necessary to determine the nodal displacements in the radial and axial direc- 
tions. The nodal displacements, u^, in the radial direction are obtained by 
the method of finite element analysis. The nodal displacements, Wj , in the 
axial direction are the same for all nodes and are obtained from the condition 
of equilibrium that the net force in the axial direction on a face perpendicu- 
lar to this direction is zero. From these nodal displacements, the principal 


strains, e r , eg, and e z in the radial, circumferential, and axial directions, 
respectively, can be calculated for each element. From these strains, the 
principal stresses o r , Og , and o z can be calculated for each element. 

In this model, it is assumed that the boron continues to expand after 
deposition. In the actual filament, such an expansion could be produced by 
the diffusion of boron from the surface into the bulk. The force driving such 
diffusion is now known; however, one might conclude that it depends upon the 
chemical vapor deposition process Itself. If a boron filament is annealed in 
an atmosphere containing oxygen, which removes boron from the surface, the 
filament contracts (Ref. 4); however, it neither contracts or elongates if an- 
nealed in an oxygen-free nitrogen atmosphere (Ref. 4) which does not remove 
boron from the surface. It is also observed that voids grow at the boron core 
interface when boron filaments are annealed in an atmosphere which removes 
boron from the surface, indicating that some of the boron has diffused from 
core-sheath interface to the surface and is not merely due to the removal of 
surface boron atoms. It appears that deposition of boron on the surface causes 
an inward diffusion of boron. 

The expansion is assumed to depend on the time after deposition of the 
boron according to 
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where t is the time after deposition of the boron in the given element, a 


is a reciprocal time constant, c 
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strains attained as at goes to infinity, and the double prime indicates 
that these are plastic strains. For large a, the plastic strains for a 
given element approach their maximum values before very many additional 
elements are added. This approaches Witucki's model in which the boron 
is deposited in a strained condition. For small a, boron in the first de- 
posited elements continues to expand even as the final elements are deposited. 

Due to these plastic deposition strains, stresses start to increase in 
the elements causing creep to take place. In this model, the following stress- 


dependent creep equations are used to calculate the creep strains for each 
element. Following Lin (Ref 5), the effective creep rate c* is given by 

e* » Bo* 0 (2) 

where B and n are constants and a*, the effective stress, is 

°* - ~~ f(°r - °e) 2 + <°e - c z> 2 + (° z - a r ) 2 l 1/2 (3) 

y/2 

The principal creep strains can then be defined in terms of e* and 0 * as 
(Ref. 5) 



where At is the time required to deposit the outermost element and is pro- 
portional to r| - Equations (4) give the principal creep strains f ,r 

an element in terras of the principal stresses for that element. 

RESULTS 

In this section, the dependence of the calculated residual stresses and 
filament elongation on the parameters in the deposition strain equations 
(Eqs . (1)) and the creep strain equations (Eqs. (2)) are shown. In these 
calculations, the elastic constants used for the carbon core and boron sheath 
are shown in Table I. 

In order to keep the next three figures simple, only the dependence of 
the axial component of the calculated residual stresses on the parameters in 
deposition strain and creep equations is shown. 

In Figure 2, the effect of varying B in the creep equations (Eq. (2)) 
on the axial residual stress, o z , is shown. For zero creep (B = 0.0) the 


axial component varies from -2.95 GPa (-A28 ksi) near the surface to 6.67 GPa 
(967 ksi) near the core. As B is increased from zero, the amount of creep 
increases with a comparable reduction in the residual stresses aB shown by the 
curves for B - 3, 6, 9, 15, and 200* lO 1 ^. 

In Figure 3, the effect of varying a in the deposition strain equations 

(Eq. (1)) on the axial residual stress, o z , is shown. Note that changes in o 

produce large changes in the shape of the curves for a z . In particular, the 
minimum just inside the outer surface changes from a broad minimum for small 
a to a narrow minimum for large o. In the limit as a goes to infinity, 
this minimum disappears (not shown in the figure). As a goes to zero, the 
surface value for c_ goes to zero. 

In Figure A, the effect of simultaneously varying the parameters z q q and 
e oz (Eq. (1)) on the axial residual stress, o z , is shown. Here e or * 0.0, and 
e 0 Q is set equal to e oz . The effects of varying this parameter are small 

near the core but become larger near the surface making large changes in the 

minimum near the surface. 

In Figure 5, the effect of varying B in the creep equations (Eq. (2)) 
on the elongation of the filament is shown. This shows the percent elongation 
versus radius of the filament during the filament growth. The largest value 
of B shows the largest reduction in elongation. 

In Figure 6, the effect of varying a in the deposition strain equation 
(Eq. (1)) on the elongation during filament growth is shown. As a increases, 
the elongation for a given radius also Increases except for the larger values 
of a in which the elongation decreases for increasing values of a for 
values of r near the largest values shown. 

In Figure 7, the effects of simultaneously varying e 0 q and e 02 (Eq. (1)) 
on the elongation during filament growth is shown. Here increasing e o6 and 
e QZ markedly Increases the elongation. 

Comparison with Experimental Results 

The parameters e Q0 , e QZ , a, B, and n were adjusted to calculate a a z 
versus r curve which closely matched the experimental curve for a 102 pm 
(A.O mil) diameter B/C filament. This is shown in Figure 8 along with an ex- 
perimental curve. The experimental curve was obtained by the method described 
in Ref. 2. Several points should be made regarding these parameters. The 


minimum in the experimental curve for o 2 near the surface is reasonably 
reproduced in the calculated curve when values of a are near 0.3. Based on 
this model) this means that the boron is not merely deposited in a strained 
state but that it continues to expand after deposition. The deep minimum at 
the core-sheath interface in the experimental curve could not be duplicated in 
the calculated curves. This experimentally observed minimum is probably the 
result of the small amount of interdiffusion between the carbon core and the 
boron sheath which was not modeled in the calculation. The lower value for 
the core residual stress in the experimental curve may result from the low 
shear strength of the pyrolytic graphite coating applied to the carbon sub- 
strate to reduce core fracture during boron deposition. Also, this was not 
considered in this model. 

In Figure 9, the calculated curves for all three components of the residual 
stresses, o r , Og, and a z are shown using the same values for the parameters 
as for the calculated curve in Figure 8. A comparison of the calculated curves 
for o r and Og with experiment cannot be made because of the lack of com- 
parable experimental data. However, the boundary condition at the surface re- 
quires that o r goes to zero as it does in this figure. 

In Figure 10, the calculated elongation results and the experimental elonga- 
tion curves of Mehalso (Ref. 1) for B/C are shown. The values of parameters 
for the calculated curve are the same as in Figures 8 and 9. The experimental 
data are for deposition temperatures of 1200° C and 1300° C are shown. Near 
the final radius, the slope of the calculated curve is decreasing more rapidly 
with radius than are the experimental curves. However, the calculated total 
elongation agrees with the experimental value. 

CONCLUSIONS 

The overall good agreement between the calculated and experimental curves 
for both the axial residual stress and elongation curves leads to the following 
conclusions : 

1. The minimum near the surface in the experimental curve for o_ is rea- 
sonably reproduced in the calculated curve for values of a near 0.3, This 
means that chemical vapor deposited boron continues to expand for some time 
after deposition. Such an expansion could result from boron diffusing from the 
surface into the bulk filling vacancies or adding additional planes of atoms. 
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2. Boron at Che temperature of deposition muBt creep substantially as 
Indicated by the large value of B required in Kquation (2). 


SUMMARY 

A three-dimensional model of the chemical vapor deposition of boron on a 
carbon substrate was developed in which an expansion of the deposited boron 
and subsequent creep of the boron were included. The method of finite element 
analysis was used to obtain values of residual stresses o r , Og , and o z as a 
function of radial location , r, in the completed filament and the elongation of 
the filament as a function of its growth radius. A wide range of values for 
parameters in the deposition strain equations and creep equations were used. 
Comparison of the calculated results with experimental results for the axial 
residual stresses and elongation of B/C filaments Indicate that the boron 
continues to expand for some time after deposition and that the boron creeps 
which limits the build-up of residual stresses. 
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